DN interaction in a color-confining chiral quark model 
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We investigate the low-energy elastic DN interaction using a quark model that confines color and 
realizes dynamical chiral symmetry breaking. The model is defined by a microscopic Hamiltonian 
inspired in the QCD Hamiltonian in Coulomb gauge. Constituent quark masses are obtained by 
solving a gap equation and baryon and meson bound-state wave functions are obtained using a 
variational method. We derive a low energy meson-nucleon potential from a quark-interchange 
mechanism whose ingredients are the quark-quark and quark-antiquark interactions and baryon 
and meson wave functions, all derived from the same microscopic Hamiltonian. The model is 
supplemented with (a, p, oj, do) single-meson exchanges to describe the long-range part of the 
interaction. Cross-sections and phase shifts are obtained by iterating the quark-interchange plus 
meson- exchange potentials in a Lippmann-Schwinger equation. Once model parameters in meson 
exchange potential are fixed to describe the low-energy experimental phase shifts of the K + N and 
K N reactions, predictions for D°N and D~ N reactions are obtained without introducing new 
parameters. 



I. INTRODUCTION 

The interaction of heavy-flavored hadrons with ordi- 
nary hadrons like nucleons and light mesons is focus of 
great contemporary interest in different contexts. One 
focus of interest is in experiments of relativistic heavy 
ion collisions. In heavy ion experiments, charm and bot- 
tom quarks are produced in the initial stages of the col- 
lision by hard scattering processes. Since they are much 
heavier than the light partons making up the bulk of the 
matter produced in the collision, it is likely that they will 
not equilibrate with the surrounding matter and there- 
fore they might be ideal probes of properties of the ex- 
panding medium - for a recent review on the subject 
and an extensive list of references see Ref. [1]. However, 
the heavy quarks will eventually hadronize and informa- 
tion on the medium will be carried out of the system by 
heavy-flavored hadrons. In their way out of the system, 
the heavy-flavored hadrons will interact with the more 
abundant light-flavored hadrons and a good understan- 
ding of the interaction is crucial for a reliable interpre- 
tation of experimental data. Another focus of interest 
is an exciting physics program that will be carried out 
with the 12 GeV upgrade of the CEBAF accelerator at 
the Jefferson Lab in the USA and the construction of 
the FAIR facility in Germany. At the Jefferson Lab 
charmed hadrons will be produced by scattering elec- 
trons off nuclei and at FAIR they will be produced by 
the annihilation of antiprotons on nuclei. One particu- 
larly exciting perspective is the possibility of creating new 
exotic nuclear bound states by nuclei capturing charmo- 
nia states like J/^S and r\ c [2-5], or heavy-light D and 
D* mesons [6-8]. In addition, the quest for experimen- 
tal signals of chiral symmetry restoration in matter is a 
subject of immense current interest and open-charm D 
and D* mesons are expected to play an important role 
in this respect - for a recent review, see Ref. [9]. In D 
mesons, properties of the light constituent quarks deter- 
mined by dynamical chiral symmetry breaking (DxSB), 



like masses and magnetic moments, are sensitive to the 
surrounding environment and changes in such properties 
will impact the structure and interactions of the mesons 
in medium. Evidently, a prerequisite for reliable predic- 
tions of modifications of hadron properties in medium 
is a good understanding of the free-space interaction of 
these hadrons. Knowledge of the interaction in free space 
is also essential for guiding future experiments aiming at 
producing exotic bound states and measuring in-medium 
hadron properties. 

The present paper is concerned with the low-energy in- 
teraction of D = (D°,D~) mesons with nucleons in free 
space. There is a complete lack of experimental infor- 
mation on this reaction, all of what is presently known 
about the interaction in free space has been apprehended 
from model calculations based on hadronic Lagrangians 
motivated by flavor SU(A) symmetry [10-12], models us- 
ing hadron and quark degrees of freedom [13] and heavy- 
quark symmetry [14, 15], and the nonrelativistic quark 
model [16]. Although the first lattice QCD studies of the 
interaction of charmed hadrons like J /iff and rj c with nu- 
cleons are starting to appear in the literature [17-19], the 
DN interaction does not seem to have been considered by 
the lattice community. In view of this situation, the use 
of models seems to be the only alternative for making the 
urgently needed predictions for e.g. low-energy cross sec- 
tions of reactions involving charmed hadrons. However, 
as advocated in the works of Refs. [13, 20-22], minimally 
reliable predictions of unknown cross-sections need to be 
founded on models constrained as much as possible by 
symmetry arguments, analogies with other similar pro- 
cesses, and the use of different degrees of freedom. In 
the specific case of the DN reaction, Ref. [13] extended 
a very successful model for the KN reaction [23], in 
which the long-range part of the interaction is described 
within a meson-exchange framework [24, 25] and the 
short-distance part is described by a quark-interchange 
mechanism from one-gluon-exchange (OGE) of a nonre- 
lativistic quark model (NRQM) [26-31]. The model of 
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Ref. [23] describes the available low-energy experimen- 
tal data for the KN reaction and was used to set limits 
on the width of the hypothetical G + (1540) pentaquark 
state [32] . For the DN reaction, the model predicts cross- 
sections that are on average of the same order of magni- 
tude but larger by a factor roughly equal to 2 than those 
of the analogous K + N and K°N reactions for centcr- 
of-mass kinetic energies up to 150 MeV. Two interesting 
findings of the study of Ref. [13] are noteworthy: (1) 
quark interchange contributes about the same amount 
as meson exchanges to the DN s-wave phase shifts, and 
(2) among the meson exchanges processes, scalar (<r) and 
vector (co,p) are the most important contributors - sin- 
gle A c , S c baryon-exchange diagrams and higher-order 
box diagrams involving D*N, DA, and D*A intermedi- 
ate states contribute very little. Recall that single-pion 
exchange is absent in this reaction. The same model was 
also used to examine the possibility to extract informa- 
tion on the DN and DN interactions in an antiproton 
annihilation on the deutcron [20] . 

The fact that quark-interchange plays a prominent 
role in the DN reaction is very significant to the quest 
of experimental signals of chiral symmetry restoration 
via changes in the interactions of D mesons in matter. 
As said, such changes would be driven at the microscopic 
level by modifications of the properties of the light u and 
d constituent quarks. However, within a NRQM there 
is no direct way to link DxSB in medium and the effec- 
tive hadron-hadron interactions, since constituent quark 
masses and microscopic interactions at the quark level 
are specified independently in the model [33]. Also, 
any temperature or density dependence on constituent 
quark masses has to be postulated in an ad hoc man- 
ner [34, 35] to account for effects of DxSB restoration. 
As a step to remedy this limitation of the NRQM, in 
the present paper we use a model that realizes DxSB in 
a way that the same microscopic interaction that drives 
D%SB also confines the quarks and antiquarks into color 
singlet hadronic states, and in addition is the source of 
the hadron-hadron interaction. The model is defined by 
a microscopic Hamiltonian inspired in the QCD Hamil- 
tonian in Coulomb gauge, in that an infrared divergent 
interaction models the full non-Abelian color Coulomb 
kernel in QCD and leads to hadron bound states that 
are color singlets - the Hamiltonian also contains an in- 
frared finite interaction to model transverse-gluon inter- 
actions and leads to, among other effects, hyperfine split- 
tings of hadrons masses [36-38] . We implement an ap- 
proximation scheme that allows to calculate with little 
computational effort variational hadron wave functions 
and effective hadron-hadron interactions that can be it- 
erated in a Lippmann-Schwinger equation to calculate 
phase shifts and cross-sections. An early calculation of 
the KN interaction within such a model, but using a 
confining interaction only was performed in Ref. [39]. 

The paper is organizes as follows. In the next Section 
we present the microscopic quark-antiquark Hamiltonian 
of the model. We discuss DxSB in the context of two 



models for the infrared divergent potentials that mimic 
full non-Abelian color Coulomb kernel in QCD, and ob- 
tain numerical solutions of the constituent quark mass 
function for different current quark masses. In Section III 
we discuss a calculation scheme for deriving effective low- 
energy hadron-hadron interactions within the context of 
the Hamiltonian of the model. A low-momentum expan- 
sion of the quark mass function is used to obtain va- 
riational meson and baryon wave functions and explicit, 
analytical expressions for the effective meson-baryon po- 
tential. In Section IV we present numerical results for 
phase shifts and cross sections for the K + N and K°N 
and D°N and D~N reactions at low energies. Initially 
we use the short-ranged quark-interchange potential de- 
rived within the model and add potentials from one- 
meson exchanges to fit experimental s— wave phase shifts 
of the K + N and K °N reactions. Next, without introduc- 
ing new parameters, we present the predictions of the 
model for the D°N and D~N reactions. Our conclu- 
sions and perspectives are presented in Section V. The 
paper includes one Appendix, that presents the meson 
Lagrangians and respective one-meson exchange poten- 
tials. 



II. MICROSCOPIC HAMILTONIAN AND THE 
CONSTITUENT QUARK MASS FUNCTION 

The Hamiltonian of the model is given as 

H = H + H lnt , (1) 

where H and H int are given in terms of a quark field 
operator $f(x) as 

H = J dx^(x)(-ia- V + pm)^{x), (2) 

and 

Hut = ~\ J dxdyp a (x)V c (\x-y\)p a (y) 

+ \ j dxdyJ?{x)D^{\x-y\)J*{y). (3) 

In the above m is the current-quark mass matrix of the 
light I = (it, d), strange s, and charm c quarks: 

/ m u \ 

m = n 7 ° ! l (4) 

m s w 
\ m c J 

and p a (x) is color charge density 

p a (x) = &(x)T a V(x), (5) 

Jf{x) is the color current density 

Jf{x) = &(x)T a ai$(x), (6) 
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with T a = A Q /2, where A a are the SU(3) Gell-Mann 
matrices. Vc and are the effective Coulomb and 
transverse-gluon interactions; the transversity of D 4 - 7 im- 
plies 



D ij {\x-y\) 



S 13 - -^r- ) D T {\x - y\). (7) 



The problem of DxSB with such an Hamiltonian has 
been discussed in the literature since long time in 
Bardeen-Cooper-Schriefer (BCS) mean-field type of ap- 
proaches, via Bogoliubov-Valatin transformations or 
Dyson-Schwinger equations in the rainbow approxima- 
tion [40-54]. For our purposes in the present paper, it 
is more convenient to follow the logic of the Bogoliubov- 
Valatin transformations. In this approach, DxSB is char- 
acterized by a momentum-dependent constituent- quark 
mass function M(k), so that the quark field operator of 
a given color and flavor can be expanded as 



/dk 
— Iu4k)q4k) + v s (k)qi(-k)}e^ 



where u s (k) and v s (k) are Dirac spinors: 



u s (k) = 



v s (k) = 



IE(k) + M(k) 
2E(k) 



(8) 



(9) 



E(k)+M(k) 



crk 



lE(k) + M(k) / ~ E{k)+M(k) 



2E{k) 



1 



(10) 



with E(k) = [k 2 + M 2 {k)] 1 / 2 , Xs is a Pauli spinor, Xs = 
—ia 2 x* S i and q\{k), qt{k) 7 q s {k), and q s (k) are creation 
and annihilation operators of constituent quarks; q s (k) 
and q s (k) annihilate the vacuum state 



q s (k)\Q)=0, 



q s (k)\Q)=0. 



(11) 



For m = 0, the Hamiltonian is chirally symmetric, but 
is not symmetric, ^ 0. 

Substituting in Eqs. (2) and (3) the expansion of "J, 
Eq. (8), and rewriting H in Wick-contracted form, one 
obtains an expression for H that can be written as a sum 
of three parts: 



H = £ + H 2 +H A , 



(12) 



where £ is the c-number vacuum energy, and H 2 and H4 
are normal-ordered operators, respectively quadratic and 
quartic in the creation and annihilation operators. The 
mass function M(k) is determined by demanding that H 2 
be diagonal in the quark operators. This leads to the gap 
equation for the constituent quark mass function Mf(k) 
of flavor /: 



M f (k) = mf + 



dq 

(2^ 



Fj 1 \k,q)V c (\k-q\) 



+ 2Gf(k iq ) D T (\k-q\) 



(13) 



where m/ is the current quark mass and 
(fe ' 9) " EM E f (q)k k 9 ' 



(14) 



Gy>(k,q) 



Ef(q) 

M f (k)q(k-q-k 2 )(k-q-q 2 ) 



+ 



E f {q) k 



kq\k — q\ 



(15) 



and Vc(k) and Dr(k) are the Fourier transforms of 
Vcd 33 !) and Dt(\x\). The quadratic Hamiltonian H 2 is 
given by 

H 2 = J2 I dks f (k) [ql f (k)q sf (k) + qt f (k)q sf (k) 

(16) 

where e / (k) is the constituent-quark single-particle ener- 
gy of flavor /, given by 



k 2 + m f M f (k) 



+ ■ 



Ef(k) 
dq 
? 



r _dq 

J (27T 



T(k,q)V c (\k-q\) 



+ 2Gf\k,q)D T (\k-q\) 



(17) 



with 



F ? , ^ = mm + wm- (18) 



M,(k) M,(q) 



k q 



and 



G^ 2) (k a) = "^ JV "> IZliyiA 4. ^" 1 Z ' v ~ 
f { > q> E f (k) E f (q) + E f (k)E f (q)\k 

The four-fermion term H A is simply the normal-ordered 
form of H int : 



Mf(k)M f (q) (k-q-k 2 )(k-q-q 2 ) 



Ha — : Hi, 



(20) 



The color confining feature of the Hamiltonian will be 
discussed in the next Section. 

To solve the gap equation on needs to specify the in- 
teractions Vc and Dt- For the confining Coulomb term 
Vc, we use two analytical forms, that we name Model 1 
and Model 2, to assess the sensitivity of results with res- 
pect to Vc- Model 1 is a parametrization of the lattice 
simulation of QCD in Coulomb gauge of Ref. [55]: 



V c {k) 



8iro~o 



Old 



fc 4 



4ttC 



(21) 



with ctcouI = (552 MeV) 2 and C = 6. Model 2 for V c was 
used in recent studies of glueballs [56] and heavy hybrid 
quarkonia [57]; it is written as 



V c (k) = V l (k) + V s (k), 



(22) 



1 „ 

where 






Model 1 : 


(23) 


1000 






with 




M(k) [MeV] 

o 
o 




a(k) = -, r-. 

/33/2 1 n 3/2 ( C +fcVA| CD ) 


(24) 




. b 


The parameters here are Aqcd = 250 MeV, Z = 
c = 40.68, and j3 — 121/12. For the transverse- 
interaction Dt we use 


5.94, 
gluon 
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\° ■ 



D T {k) 



{k 2 + m 2 )ln LA2 (T + k 2 /m 2 ) 



(25) 



This choice is guided by previous studies of spin-hypcrfinc 
splittings of meson masses [38] using an Hamiltonian as 
in the present work. Moreover, the Yukawa term multi- 
plying the log term is used to conform lattice results that 
the gluon propagator in Coulomb gauge is finite in the 
infrared [55, 58]. Further ahead we will discuss the im- 
pact of this choice of infrared behavior on our numerical 
results. Parameters here are m g = 550 MeV, m = m g /2, 
t = 1.05, and olt — 0.5. We use the same Dt for both 
models. We note that we could have used a log running 
similar to the one in Vc, Eq. (24), but results would not 
change in any significant way. 

We have solved the gap equation in Eq. (13) by ite- 
ration. The angular integrals can be performed analyti- 
cally, but special care must be taken with the strongly- 
peaked confining term l/|fe — q| 4 at q s» k in the nu- 
merical integration over q. There is no actual divergence 
here: the terms M(q)/E(q) and M(k)/E(q) (q/k) k ■ q 
in Eq. (14) cancel exactly when q = k, but this cancel- 
lation can be problematic in the numerical integration. 
This problem can be handled in different manners, like 
using a momentum mesh containing node and half-node 
points with e.g. k at nodes and q at half-nodes [44], or in- 
troducing a mass parameter /im such k 4 — » (k 2 +/i 2 ) 2 in 
Eq. (21) and varying pm until results become indepen- 
dent of pin [59] . In the present paper we use a different 
method [60]; we add a convenient zero to the gap equa- 
tion: 



dq 



1 



M(k) 1 
(2tt) 3 \k-q\ 4 E(k) k 



k- (k — q) 



(26) 



and rewrite F^'(k,q) in Eq. (13) as 

FW(fe,g)^FW(fe,g)-|Hlfc.( fe _ g) 

= (M(q) M(k)\ 
\E(q) E(q)J 



fM(k) M(k)\ 



(27) 



so that angle-independent and angle-dependent terms 
vanish independently when q = k. This feature makes 
the numerical cancellation of the divergences very stable. 
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FIG. 1. Constituent quark mass M(k) as a function of the 
momentum for different values of current-quark masses. 



In Fig. 1 we present solutions M(k) of the gap equa- 
tion for different values of current-quark masses. The 
dramatic effect of D^SB mass-generation is seen clearly 
in the figure: for the fictitious limit of a zero current- 
quark mass (solid line), the mass function M(k) acquires 
a sizable value in the infrared, M(k = 0) sa 250 MeV. 
In the ultraviolet, the mass function runs logarithmi- 
cally with A:, in a manner dictated by the running of the 
microscopic interactions. Since the model interactions 
used here have a different ultraviolet running from the 
one dictated by perturbative QCD, the logarithm run- 
ning of the quark-mass function must be different; the 
interactions here fall off faster than those in QCD. One 
consequence of this is that the momentum integral over 
the trace of the quark propagator, which gives the quark 
condensate, does not run and is ultraviolet finite in the 
present case; we obtain (-('i"J')) 1 / 3 = 280 MeV. For 
nonzero current-quark masses, Fig. 1 shows that the ef- 
fect of mass generation diminishes as the value of m/ 
increases - here m u = ma = 10 MeV, m s — 150 MeV 
in both models, and m c = 950 MeV in Model 1 and 
m c = 600 MeV in Model 2. On the logarithm scale, one 
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sees that the mass function varies substantially in the 
ultraviolet, from k ~ 3Aqcd = 750 MeV onwards. It is 
worth noticing that although the momentum dependence 
of the mass function charm quark is much less dramatic 
than for light quarks, there is still a significant dressing 
effect, M c (0)/to c <~ 1.5 — 3, similarly to what is found 
with covariant Dyson-Schwinger equations [61-63]. 

Finally, we comment on the impact of the momentum 
dependence of Dx{k) in the infrared on the mass func- 
tion. In the complete absence of D T , using only Vc with 
the parameters above, M(0) would be of the order of 
100 MeV in the chiral limit. On the other hand, using a 
D T {k) that vanishes at k = 0, like in the form of a Gribov 
formula [64], would give M(0) ~ 200 MeV in the chiral 
limit. This is in line with the recent finding of Ref. [65], 
which studies D^SB in a framework based on a quark 
wave functional determined by the variational principle 
using an ansatz which goes beyond the BCS-type of wave 
functional. 



III. BARYON-MESON INTERACTION 

In this Section we set up a calculation scheme for 
deriving effective low-energy hadron-hadron interactions 
with the Hamiltonian discussed above. We seek a scheme 
that makes contact with traditional quark-model calcula- 
tions and can be systematically improved with some com- 
putational effort. Since the early 80's, the great majority 
of calculations of hadron-hadron scattering observables in 
the quark-model are based on methods to handle cluster 
dynamics adapted from nuclear or atomic physics - for 
reviews on methodology and fairly complete lists of refe- 
rences see Refs. [66-68]. Such methods require a micros- 
copic Hamiltonian and hadron bound state wave func- 
tions given in terms of quark degrees of freedom. Within 
the context of the model discussed in the previous Sec- 
tion, the starting point is the microscopic Hamiltonian of 
Eq. (1), with quark field operators given in terms of the 
constituent quark mass function M(k) obtained from the 
numerical solution of the gap equation in Eq. (13). 

The sector of the Hamiltonian relevant for the elastic 
meson-baryon interaction can be written in a compact 
notation as 

Hqq = e(M) qIQ/m + qUv + ^V qq (pv; op) q\q\q p q a 

+\ v qq{w, G P) qIqIqpQo + Vqqinv; op) q\qtq p q a , (28) 

where the indices p, v, p, and a denote collectively the 
quantum numbers (orbital, color, spin, flavor) of quarks 
and antiquarks. The first two terms in Eq. (28) are the 
quark and antiquark single-particle self-energies coming 
from Eq. (16), and V qq , V qq , and V qq are respectively the 
quark-quark, quark- antiquark and antiquark- antiquark 
interactions from H4 in Eq. (20). The one-baryon state in 
the BCS approximation can be written in the schematic 



notation as [47, 51, 52] 

\a) = Bl\Q) = -±= V£™ qIqUUM, (29) 

with ^£iM2M being the Fock-space amplitude, with a de- 
noting the (orbital, spin, flavor) quantum numbers of the 
baryon. Likewise, the one-meson states is written as 

|o) - Mt|fi) = Cg+gJlfi), (30) 

where <f>% v is the corresponding Fock-space amplitude, 
with a representing the quantum numbers of the meson. 
The Fock-space amplitudes tp and (f> can be obtained by 
solving a Salpeter-type of equation [47, 51, 52]. 

Given a microscopic Hamiltonian H qq and hadronic 
states I a) as in Eqs. (28)-(30), an effective low-energy 
meson-baryon Vjfs(a6, cd) potential for the process 
Meson(a) + Baryon(6) — > Meson(c) + Baryon(d) can be 
written as [67] 

V MB (ab,cd) =-Z4>T^7^V qq ^o P )4>TV h ^ 

- 3 ^r^T^v qq {pv- Mc^r" 3 

- 6 ^i,* d V ^V qq {pv; ap)4>r^r^ 

-60; w ^>;^)€T' w . (si) 

where </>^, ip^ va , • • • are the meson and baryon Fock- 
space amplitudes of Eqs. (29) and (30), and V qq , V qq , and 
V qq are respectively the quark-quark, quark-antiquark 
and antiquark-antiquark interactions in H qq of Eq. (28). 
The expression in Eq. (31) is completely general, in that 
it is valid for any low-energy meson-baryon process for 
which the baryon and meson Fock-space states and H qq 
are as given above. It can be iterated in a Lippmann- 
Schwinger equation to obtain scattering phase shifts and 
cross sections - for details, see Ref. [67]. There is, how- 
ever, one difficulty: Vmb involves multidimensional in- 
tegrals over internal quark and antiquark momenta of 
products of baryon and meson wave functions and prod- 
ucts of Dirac spinors u(k) and v(k) that depend on the 
quark mass function M(k). Although this is not a major 
difficulty, an approximation can be made noting that the 
bound-state amplitudes tp and 4> are expected to fall off 
very fast in momentum space for momenta larger than 
the inverse size of the hadron, and therefore only low- 
momentum quark and antiquark processes contribute in 
the multidimensional integrals. In view of this, and the 
fact that M{k) changes considerably only at large mo- 
menta, a natural approximation scheme to simplify the 
calculations without sacrificing the low-energy content of 
the effective interaction, is to retain the first few terms 
in the low momentum expansion for the mass function 
Af(fc): 

M(k) = M + M'(O) k + ^M"(0) k 2 + ■ ■ ■ , (32) 

where M'(fc) = dM(k)/dk and M"(k) = d 2 M(k)/dk 2 . 
In particular, retaining terms up to 0(k 2 /M 2 ) in the ex- 
pansion, it is not difficult to show that the Dirac spinors 
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u s (k) and v s (k) in Eqs. (9) and (10) become 

1 - 



u s (k) 



v s (k) 



Xs ) 




(33) 



(34) 



Using these spinors in Eq. (3), the expressions one ob- 
tains for V qq , V q q, and Vqq are very similar to those of 
the Fermi-Breit expansion of the OGE interaction [69]. 
There is, however, one important difference here: while 
for the OGE one has V c (k) = D T (k) « 1/fc 2 , in the 
present case Vc(k) and Dx{k) arc different and repre- 
sent very different physics; Vc is a confining interaction 
and Dt is a (static) transverse gluon interaction. 

The evaluation of multidimensional integrals can 
be further simplified using the variational method of 
Refs. [47, 51] with Gaussian ansatze for the bound- 
state amplitudes ip and <f>, instead of solving numerically 
Salpeter-type of equations. Specifically: 



ipp(ki,k 2 , k 3 ) = 8(P - fei - k 2 - k 3 ) 



3/4 



e -EL 1 (*,-P/3)/2« 2 ](35) 



>(k 1 ,k 2 ) = 5(P-k 1 -k 2 ) 



7T/3 2 



3/4 



-(Mife 1 -M 2 fe 2 ) 2 /8/3 2 



(36) 



where a and j3 are variational parameters, P is the 
ccntcr-of-mass momentum of the hadrons and 



Mi 



2M S 



(M q + Mg) ' 



M 2 



2M„ 



(M q + Mg 



(37) 



with M q and Mq being the zero-momentum constituent- 
quark masses. The variational parameters are deter- 
mined by minimizing the hadron masses: 



M„ 



(a\{H 2 + H 4 )\a'} 



(a\a'} 



(38) 



where we left out the constant £, defined in Eq. (12), 
which cancels in the hadron mass differences - see Table I. 

In this work we consider the elastic scattering of the 
pseudoscalar mesons (K + ,K°) and (D°,D~) off nu- 
cleons, with both the mesons and nucleons in their 
ground states. Using the Gaussian forms for the am- 
plitudes tp and <f>, Eqs. (35) and (36), one obtains for the 
nucleon mass, to at, and for the pseudoscalar meson mass, 
mp, the following expressions: 



Mm = 3 



+ 



+ 



3/2 



dk 



k 2 MAk) 
+ mi 



Ei(k) 1 E t (k) 



( 3 \ 3/2 f dkdg 
\2na 2 ) J (2tt) 3 



/' dkdq 
(2^ 



2F^(k,q) + 3C N e-^^ 2a2 



Vc{\k-q\)e-^/^ 2 



2na 2 



3/2 . 



4Gf (fc,<?) 



(fc ~ <?) 2 e -{k- q flio? 
3M, 2 



D T (\k-q\)e 



-k 2 /c 



(39) 



M P = 



7T/3 2 



3/2 



dk 



7T/3 2 



3/2 



Et(k) 
dkdq j 2 
(2^)3 I 3 



mi 



Mi(k) k 2 



m h 



M h (k) 



E h {k) 



-fc 2 //3 2 



F l {2 \k,q)+Fl i 2 \k,q) +C P e-^/ 2 ^ }V c (\k-q\)e 



-fe 2 //3 2 



+ ■ 



J_\ 3/2 fdkdq\ G (2) ( a) + G (2) (k Q) . l (fc-g) 2 r -(k- Q f/ 20 * 
nP 2 ) J (2tt)3 G; ^1)+ G h [k,q)+ 2 e 



D T (\k-q\)e 



-k 2 /p 2 



(40) 



where h and G\l(k,q) are given in Eqs. (18) 

and (19) and the indices I and h refer to light and heavy 
flavors, I = (u,d) and h — (s,c). The matrix ele- 
ments of the color matrices, Cat = (N\T a T a \N) and 
C P = (P\T a (-T a ) T \P), in the terms proportional to 
the Coulomb potential Vc - that come from H4 - are 



written in a way to emphasize the color-confinement fea- 
ture of the model [47]: for k = q, these terms are diver- 
gent, unless the color matrix elements are such that the 

( 2) 

corresponding expressions vanish. The terms F^(k,q) 
come from H 2 - which is diagonal in color. For baryon 
and meson color singlet states, the contributions from H 2 
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and H4 cancel exactly for k = q, since: 



Cm — 



r CiC 2 C3 e C 1 C 2 C3 

3! 



/rpa\cic 1 ^jia^c 2 c 2 



gc lC2 fic lC2 



and 



limF, (2 ft ) (fe,q) = l. 



2 

~3' 

4 
3' 



(41) 



(42) 



(43) 



Such a cancellation of infrared divergences also plays an 
important role in context of a conjectured [70] new high- 
density phase of matter composed of confined but chirally 
symmetric hadrons [71, 72]. 

Next, we obtain an explicit expression for the effec- 
tive meson- nucleon interaction Vmb, given generically 
in Eq. (31). This effective interaction is generated by 
a quark-interchange mechanism. The use of Gaussian 
wave functions is very helpful for getting a Vmb in closed 
form [67] : it can be written as a sum of four contributions 
(see Eqs. (9) and (10) of Ref. [13]), each contribution 
corresponding to a quark-interchange diagram shown in 
Fig. 2: 



1 4 

Vmb (p, p') = 2 E ^ t ^ (P- + V > & > P) 1 

i=l 



(44) 



where p and p' are the initial and final center-of-mass 
momenta, and the Vi(p,p') are given by 



Vi{p,p') = 



35 



(3 + 2g)-Ka 2 
dq 



(27T) 



v(q)e 



- aip 2 —bip' 2 +Cip-p' 



-dig +e.i-q 



(45) 



where the m, bi 7 a, di, and a involve the hadron sizes a 
and /3 and the constituent quark masses M u , M s and M c , 
and g = (a//3) 2 . Since we are using the same Fock-space 
amplitudes tp and cf) as those used in Ref. [13], the expres- 
sions for at, hi, ■ ■ ■ are the same as there. It is, howe- 
ver, important to note that the essential difference here 
as compared to Ref. [13] is that the constituent quark 
masses M u , M s , and M c , the width parameters a and 
/?, and the mcson-baryon interaction are all derived from 
the same microscopic quark-gluon Hamiltonian. Another 
very important difference here is v(q): while in Ref. [13] 
it comes from the OGE, here v(q) = Vc(q) for the spin- 
independent interaction and v(q) = 2</ 2 /(3M;M') Dx(q), 
for the spin-spin interaction, with M' = Mi for dia- 
grams (a) and (c), and M' = Mh for diagrams (b) and 
(d). The LOi are coefficients that come from the sum over 
quark color-spin-flavor indices and combinatorial factors, 
whose values for the DN interaction are given in Table 3 
of Ref. [13]. The corresponding coefficients for the KN 
reaction can be extracted from that table by identifying 
D° with K+ and D~ with K°. Note that V M b{p,p') is 



symmetric under interchange of p and p' and therefore 
possess post-prior symmetry [67, 73], a feature that is 
not always satisfied when using composite wave functions 
that art not exact eigenstates of the microscopic Hamil- 
tonian [74]. As mentioned earlier, the quark- interchange 




(a) 



(b) 



(c) 



(d) 



FIG. 2. Pictorial representation of the quark-interchange pro- 
cesses that contribute to a meson-baryon interaction. The 
curly lines represent the interactions Vc and Dt- 



mechanism leads to an effective meson-baryon potential 
that is of very short range. It depends on the overlap 
of the hadron wave- functions and contributes mostly to 
s- waves. As shown in Ref. [23], in order to describe 
experimental data for the K + N and K°N reactions, 
light-meson exchange processes are required to account 
for medium- and long-ranged components of the force. 
Quark-interchange generated by OGE accounts only to 
roughly 50 % to the experimental s-wave phase shifts 
and vector (u> and p) and scalar exchanges are crucial 
for the correct description of this wave as well as higher 
partial waves [23]. Moreover, in order to describe the cor- 
rect isospin dependence of the data, it is essential to in- 
clude the exchange of the scalar- isovector clq meson [23] . 
With these facts in mind, we follow Ref. [13] and include 
meson-exchanges in the DN system. As in that refe- 
rence, we parameterize correlated 7T7t exchange in terms 
of a single c-meson exchange - this is not a bad appro- 
ximation for the 1 = channel, but for I = 1 channel 
it underestimates the total strength by 50% [13]. In Ap- 
pendix A we present the effective Lagrangians densities 
and corresponding one-meson-exchange amplitudes. 



IV. NUMERICAL RESULTS FOR PHASE 
SHIFTS AND CROSS SECTIONS 

As mentioned previously, not much is known experi- 
mentally about the DN interaction at low energies. In 
view of this and in order to have a comparison standard, 
we also consider within the same model, without chan- 
ging any parameters besides the current quark masses, 
the K + N and K°N clastic processes, for which there 
are experimental data. It is important to reiterate that 
the effective short-range meson-baryon potential derived 
from quark-interchange driven by the microscopic inte- 
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ractions of Model 1 or Model 2 are very constrained, 
in the sense that they are determined by the micros- 
copic interactions Vc and Dt via a chain of interme- 
diate results for hadron properties driven by D^SB. 
That is, the microscopic interactions Vc and Dt deter- 
mine the constituent quark masses Mi = (M u ,Md) and 
M/j = (M S ,M C ), and these quark masses together with 
Vc and Dt determine the hadron wave functions. The 
effective meson-baryon interaction is determined by all 
these ingredients simultaneously, since it depends expli- 
citly on the hadron wave functions, the quark masses 
and Vc and Dt- Such an interdependent chain of results 
determining the effective meson-baryon interaction is ab- 
sent in the NRQM with OGE, where quarks masses are 
independent of the microscopic quark-antiquark Hamil- 
tonian. 

Having determined the constituent quarks masses 
Mi = (M u ,M d ) and M h = (M S ,M C ), the next step is 
the determination of the variational parameters a and /3 
of the and <f> wave functions. Minimization of Afjv and 
Mp with respect to a and /3 leads to the results shown 
in Table I; also shown are the parameters used with the 
OGE model in Refs. [13, 23]. 

As expected, because the charm quark is heavier than 
the strange quark, the charmed D mesons are smaller 
objects than the kaons - recall that the meson root mean 
square radii are inversely proportional to /3. One also 
sees that hadron sizes of Models 1 and 2 are smaller than 
those used in the OGE model. Although the sizes of the 
wave functions have an influence on the degree of overlap 
of the colliding hadrons, the microscopic interactions Vc 
and Dt also play a role in the effective meson-baryon 
potentials. As we shall discuss shortly ahead, there is 
an interesting interplay between these two effects. For 
a recent discussion on the influence of hadron sizes on 
the quark- interchange mechanism, see Ref [75]. Regard- 
ing the hadron masses, one sees a discrepancy of 25% 
(Model 1) and 50% (Model 2) on the kaon-nucleon mass 
difference; this does not come as a surprise in view of 
the pseudo-Goldstone boson nature of the kaon, in that 
a BCS variational form for the kaon wave function is a 
poor substitute for the full Salpeter amplitude [5f, 52]. 
On the other hand, the DN mass difference is within 10% 
of the experimental value in both Model 1 and Model 2. 



from the quark-interchange mechanism and one-meson 
exchanges, using the method discussed in Section 2.4 
of Ref. [76]. The importance of going beyond quark- 
Born diagrams by iterating the quark-interchange po- 
tentials in a scattering equation has been stressed pre- 
viously [13, 23, 68]. For the specific case of the OGE 
quark-interchange in DN scattering, unitarization of 
the scattering amplitude by iteration in a Lippmann- 
Schwinger equation leads to a decrease of the cross- 
section at low energies by a factor of three as compared 
to the nonunitarized quark-Born amplitude [13]. 
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E cm " m N " m K t MeV ] 
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75 100 125 150 



E cm " m N " m K t MeV l 



TABLE I. Variational size parameters of the hadron ampli- 
tudes and hadron mass differences. All values are in MeV. 
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Pd 




AM DN 
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Experiment 








(443) 


(928) 


(1371) 


Model 1 


568 


425 


508 


350 


990 


1345 


Model 2 


484 


364 


423 


205 


1010 


1220 


OGE Refs. [13, 23] 


400 


350 


383.5 









We are now in position to discuss numerical results 
for scattering observables. We solve numerically the 
Lippmann-Schwinger equation for the potentials derived 



FIG. 3. KN s-wave phase shifts for isospin I — (upper 
panel) and I — 1 (lower panel) channels from the quark- 
interchange meson-baryon potential driven by the microscopic 
interactions of Model 1 (solid lines) and Model 2 (dashed 
lines). The curve with diamond symbols are the results from 
the phase shift analysis from the George Washington Data 
Analysis Center (DAC), Ref. [77]. 

In Fig. 3 we present results for the s-wave phase shifts 
for the isospin 1 = and J = 1 of the K + N and 
K°N reactions. Here we show the results derived from 
the baryon-meson potential obtained from the quark- 




FIG. 4. Same as in Fig. 3, but with meson-exchange potentials 
added to the quark-interchange meson-baryon potential. 



interchange mechanism driven by the microscopic inter- 
actions of Model 1 and Model 2. Also shown in the fig- 
ure are results from phase shift analysis from the George 
Washington Data Analysis Center (DAC) [77]. As in 
the case of the model with OGE [23], both Model 1 and 
Model 2 reproduce the experimental fact that the s— wave 
phase shifts for I = are much smaller than for 7 = 1. 
This is due to the combined effects that for 1 = 0, the 
confining interaction Vc does not contribute at all and 
Dt contributes only via diagrams c and d of Fig. 2. 

Also seen in Fig. 3 is the fact that both Model 1 and 
Model 2 provide large repulsion - the confining interac- 
tion Vc of Model 2 provides a little more repulsion than 
the Vc from Model 1. As in the case of OGE, meson- 
exchanges can be added to obtain a fair description of 
the data, as we discuss next. 

In Fig. 4 we present results for the phase shifts when 
(a, uj, p, and ao) one-meson exchanges are added to the 
quark-interchange potential. The input parameters for 
the meson-exchange potential are cutoff masses in form 
factors and coupling constants - we take the values used 



FIG. 5. DN s-wave phase shifts for 1 = (upper panel) 
and I — 1 (lower panel) channels from the quark-interchange 
meson-baryon potential driven by the microscopic interac- 
tions of Model 1 (solid lines) and Model 2 (dashed lines). 



in Refs. [13, 23], with the exception of the a and ao cou- 
plings, the product g a MM9aBB is increased by four and 
the product g aQ MM9a BB is increased by three for a rea- 
sonable description of the phases - this is because the 
quark-interchange potential gives strong repulsion. Once 
the K+N and K°N phase shifts are fitted, we use the 
same values of the cutoff masses and couplings to make 
predictions for the DN system, which we discuss next. 

In Fig. 5 we present the s— wave phase shifts for isospin 
1 = and 1=1 states of the D°N and D~N reactions. 
Results are obtained with a meson-baryon potential from 
quark-interchange driven by the interactions of Model 1 
and Model 2. Like in the similar KN system, the phases 
for the 1=1 channel are much bigger than those for 
the 1 = channel. Also, one sees that Model 2 gives a 
stronger repulsion. 

Adding (er, w, p, and ao) one-meson exchanges to the 
quark-interchange potential leads to the results shown 
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FIG. 6. Same as in Fig. 5, but with meson-exchange potentials 
added to the quark-interchange meson-baryon potential. 



in Fig. 6. Parameters of the meson-exchange potentials 
are the same used for the KN. The predictions for the 
s— wave phase shifts for DN system are qualitatively si- 
milar to the results for the KN system, but are roughly 
a factor of two larger than the latter ones. 

The quark-interchange meson-baryon potential leads 
to much smaller phase shifts for the higher partial waves. 
Like in the KN system, meson exchanges play a much 
more important role in these waves. Although not shown 
here, many important features of the meson-exchanges 
discussed in Ref. [13], as the interference pattern between 
p and u> contributions - destructive for 1 = and con- 
structive for 1=1- are seen in the present model as 
well. We do not show these detailed results here, as they 
concern the meson exchange part of the full interaction, 
but present in Fig. 7 the final predictions for DN cross 
sections. 

The results from both confining models are qualita- 
tively similar, although the results from Model 1 for the 
1 = 1 state are on average larger by a factor roughly 
equal to two. It is important to note that in the present 



FIG. 7. Predictions of DN cross sections for 1 = (upper 
panel) and 1 = 1 (lower panel) using the full meson-baryon 
potential with quark-interchange and meson exchanges. 



paper we are approximating the correlated two-pion ex- 
change contribution by single cr-meson exchange. As seen 
in Ref. [13], this seems to be a reasonable approximation 
for the 1 = channel, but underestimates the 1 = 1 
cross section by a factor roughly equal to two. We ex- 
pect a similar feature in the present case. As a final re- 
mark, we note that the results for the cross sections are of 
the same order of magnitude as those from OGE quark- 
interchange, but the shapes of the curves for the 1=1 
cross sections from Model 1 and Model 2 arc different at 
very low energies from those from OGE. This last feature 
can possibly be attributed to the drastically different mo- 
mentum dependences of Vc(k) and Vr(k) from the 1/k 2 
dependence of the OGE. 
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V. CONCLUSIONS AND PERSPECTIVES 

There is great contemporary interest in studying the 
interaction of charmed hadrons with ordinary hadrons 
and nuclei. Of particular interest is the study of D 
mesons in medium, mainly in connection with the pos- 
sibility of creating new exotic nuclear bound states by 
nuclei capturing charmonia states like J/\& and rj c [2-5], 
or D and D* mesons [6-8]. One major difficulty here 
is the complete lack of experimental information on the 
free-space interactions that would be of great help to 
guide model building - this is not the case for the si- 
milar problems involving strange hadrons. In a situation 
of lack of experimental information, one way to proceed 
in model building is to use symmetry constraints, analo- 
gies with other similar processes, and the use of different 
degrees of freedom [13, 20-22]. With such a motivation, 
we have implemented a calculation scheme for deriving 
effective low-energy hadron-hadron interactions based on 
a phenomenological Hamiltonian inspired in the QCD 
Hamiltonian in Coulomb gauge. The model Hamilto- 
nian, defined in Eqs. (l)-(7), confines color and realizes 
D%SB. The scheme makes contact with traditional cons- 
tituent quark models [33], but it goes beyond such mo- 
dels: the constituent quark masses are derived from the 
very same interactions that lead to hadron bound states 
and hadron-hadron interactions, while in the traditional 
quark models, quark masses and quark-quark interac- 
tions are specified independently. Moreover, the model 
confines color, in that only color-singlet hadron states are 
of finite energy. These features of the model are essential 
for studies seeking signals of in-medium modifications of 
hadron properties. 

The model Hamiltonian requires as input a Coulomb- 
like term Vc and a transverse-gluon interaction Dt- 
We used two analytical forms for Vc- Model 1 is a 
parametrization of the lattice simulation of QCD in 
Coulomb gauge of Ref. [55] and Model 2 was used in 
recent studies of glueballs [56] and heavy hybrid quarko- 
nia [57]. For Dt we were guided by previous studies 
of spin-hyperfine splittings of meson masses [38]. Ini- 
tially we obtained the constituent quark mass function 
Mf(k) for different flavors /. Next, we derived an effec- 
tive low-energy meson-baryon interaction from a quark- 
interchange mechanism, whose input is the microscopic 
Hamiltonian and hadron bound state wave functions 
given in terms of quark degrees of freedom [67, 68]. 
We derived explicit analytical expressions for the effec- 
tive meson-baryon interaction by using a low-momentum 
expansion of the constituent quark mass function and 
variational hadron wave functions [47, 51]. Initially we 
used the short-ranged quark-interchange potential de- 
rived within the model and added one-meson-exchangc 
potentials to fit experimental s— wave phase shifts of the 
K + N and K°N reactions. Next, without changing any 
parameters besides the strange to charm current quark 
masses, m s — > m c , we presented the predictions of the 
model for the D°N and D~N reactions. The results for 



the cross sections obtained with the present model are of 
the same order of magnitude as those from OGE quark- 
interchange, of the order of 5 mb for the 1 = 1 state and 
10 mb for the 1 = state, on average. However, the 
shapes of the curves for the 1=1 cross sections from 
Model 1 and Model 2 are different at very low energies 
from those from OGE - a feature we attribute to the dras- 
tically different momentum dependences of Vc(k) and 
Vr(k) from the 1/k 2 dependence of the OGE. 

The model can be improved in several directions. First 
of all, without sacrificing the ability of obtaining analyti- 
cal expressions for the effective meson-baryon interaction, 
one can expand the the amplitudes ip and <j> in a basis 
of several Gaussians and diagonalizc the resulting Hamil- 
tonian matrix together with a variational determination 
of the size parameters of the Gaussians. Another im- 
provement, that certainly will be necessary for studying 
in-medium chiral symmetry restoration, is to use the full 
mass function Mf (k) for the light u and d flavors, instead 
of the low- momentum expansion of Eq. (32). This is be- 
cause at finite baryon density and/or temperature the 
mass function Mf(k) for / = (u,d) will loses strength in 
the infrared [45, 80-83] and the expansion in powers of 
k 2 /M 2 evidently loses validity. However, this will have a 
slight impact on the numerics, since in this case multidi- 
mensional integrals for the determination of the hadron 
sizes and the effective meson-baryon interaction need to 
be performed numerically. The low-momentum expan- 
sion must also be abandoned in the calculation of hadron 
wave functions when the quark mass function Mf(k) in 
the infrared is much smaller than the average momentum 
of the quark in the hadron. Such a situation can happen 
when Dxik) is suppressed in the infrared, as in the study 
of Ref. [65]. 

Finally, another subject that needs careful scrutiny 
is the use of SU(4) flavor symmetry Eq. (A10), to 
fix the coupling constants in the effective meson La- 
grangians. A recent estimation [84] , using a framework in 
which quark propagators and hadron amplitudes are con- 
strained by Dyson-Schwinger equation studies in QCD, 
finds that while SU(3)-flavor symmetry is accurate to 
20%, SU(4) relations underestimate the DDp coupling 
by a factor of 5. On the other hand, a study employing a 
3 Po pair-creation model with nonrelativistic quark model 
hadron wave functions finds smaller SU(4) breakings [85]. 



Appendix A: Meson-exchange contributions 

The meson-exchange contributions are obtained from 
the following effective Lagrangian densities [78, 79]: 

Cnns(x) = g N NS^ N) (x)T ■ ^ s \x)^ N \x), (Al) 
Cnnv(x) = g NNV U {n \x)^^ n \x)t ■ ^{x) 



+ 



(~2.\ 



$ {N) {x)a tlv ^ (N) (x)t ■ d v <f><p (x)] , (A2) 
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Cpps{x) = 9pps^ p \x)^ s \x)^ p \x), 



(A3) 



Cppv(x) = igppv ^ P) (x) (d v <p( p \x)} r 

- (d v ^ p \x)) t^ p \x)] ■ 4><p(x). (A4) 

In these, i/jW denotes the nucleon doublet, <f>^ p \x) the 

charmed and strange meson doublet, cj)^\x) the iso- 
triplet of p mesons, and r are the Pauli matrices. The 
Lagrangians for the o and w mesons are obtained taking 
t — y 1 in the expressions above and in addition Ky = 
for the case of co. 

The tree-level potentials derived from the above La- 
grangian densities lead to the following expressions for 
the vector- meson exchanges (v — p,tv): 



V v {p, P ') 



9nnv gppv 



A v (ps,p's') + 



2m jv 



B u (ps,p's') 



, (A5) 



and for the scalar-meson exchanges (S = a,ao) 

As(q)u(p',s')u(p,s), 
(A6) 



V Sf i ) = 9nns 9pps a _/„\.-.r^ „/^ 

Kir* ! f ) (2 ^ 3 ^____ 



where is the nucleon mass and = (q 2 + m 2 ) 1//2 
with m the meson masses, A^(g) and As(q) are the 
vector-meson and scalar-meson propagators, u(p, s) are 
the Dirac spinors of nucleons (same expression as in 
Eq. (9), with M{k) replaced by the nucleon mass mjv), 
and the quantities and i? M are given by 



Ap(ps,p's') = u(p',s')7m u (p. s ): 



(A7) 



B^ips.p s ) = u(p , s ) ia^ q v u(p, s). (A8) 

To avoid divergences in the Lippmann-Schwinger equa- 
tion, the meson-exchange potentials are regularized phe- 
nomenologically by monopolc form factors at each vertex: 



(A9) 



where q = p' —p, mi is the mass of the exchanged meson 
and Aj is a cutoff mass. The coupling constants are fixed 
by SU(4) symmetry as in Rcf. [13]: 



9dd p - 9dDu - 9kk p = 9kklj — 



(A10) 



[16] 



R. Rapp, D. Blaschke, and P. Crochet, Prog. Part. Nucl. [17 
Phys. 65, 209 (2010) 

S. J. Brodsky, I. A. Schmidt, and G. F. de Teramond, [18 
Phys. Rev. Lett. 64, 1011 (1990). 

S. H. Lee and C. M. Ko, Phys. Rev. C 67, 038202 (2003). [19 
G. Krein, A. W. Thomas and K. Tsushima, Phys. Lett. 
B 697, 136 (2011). [20 
K. Tsushima, D. H. Lu, G. Krein and A. W. Thomas, 
Phys. Rev. C 83, 065208 (2011). [21 
K. Tsushima, D. H. Lu, A. W. Thomas, K. Saito, and 
R. H. Landau, Phys. Rev. C 59, 2824 (1999). [22 
C. Garcia-Recio, J. Nieves and L. Tolos, Phys. Lett. B 
690, 369 (2010). [23 
C. Garcia-Recio, J. Nieves, L. L. Salcedo and L. Tolos, 
Phys. Rev. C 85, 025203 (2012). [24 
R. S. Hayano and T.Hatsuda, Rev. Mod. Phys. 82, 2949 
(2010). [25 
T. Mizutani and A. Ramos, Phys. Rev. C 74, 065201 
(2006) [26 
Z. -w. Lin, C. M. Ko and B. Zhang, Phys. Rev. C 61, 
024904 (2000) . [27 
J. Hofmann and M. F. M. Lutz, Nucl. Phys. A 763, 90 
(2005) [28 
J. Haidenbauer, G. Krein, U.-G. Meifiner, and A. Sibirt- 
sev, Eur. Phys. J. A 33, 107 (2007). [29 
S. Yasui and K. Sudoh, Phys. Rev. D 80, 034008 (2009). 
C. Garcia-Recio, V. K. Magas, T. Mizutani, J. Nieves, [30 
A. Ramos, L. L. Salcedo and L. Tolos, Phys. Rev. D 79, 
054004 (2009). [31 
T. F. Carames and A. Valcarce, Phys. Rev. D 85, 094017 
(2012). [32 



K. Yokokawa, S. Sasaki, T. Hatsuda and A. Hayashigaki, 
Phys. Rev. D 74, 034504 (2006) . 

T. Kawanai and S. Sasaki, Phys. Rev. D 82, 091501 
(2010) . 

L. Liu, H. -W. Lin and K. Orginos, PoS LATTICE 2008, 
112 (2008). 

J. Haidenbauer, G. Krein, U. -G. Meissner and A. Sibirt- 
sev, Eur. Phys. J. A 37, 55 (2008). 

J. Haidenbauer and G. Krein, Phys. Lett. B 687, 314 
(2010). 

J. Haidenbauer, G. Krein, U. -G. Meissner and L. Tolos, 
Eur. Phys. J. A 47, 18 (2011). 

D. Hadjimichef, J. Haidenbauer and G. Krein, Phys. Rev. 
C 66, 055214 (2002). 

R. Biittgen, K. Holinde, A. Miiller-Groeling, J. Speth, 

and P. Wyborny, Nucl. Phys. A 506, 586 (1990). 

M. Hoffmann, J.W. Durso, K. Holinde, B.C. Pearce, and 

J. Speth, Nucl. Phys. A 593, 341 (1995). 

T. Barnes and E. S. Swanson, Phys. Rev. C 49, 1166 

(1994). 

B. Silvestre-Brac, J. Leandri and J. Labarsouque, Nucl. 
Phys. A 589, 585 (1995). 

B. Silvestre-Brac, J. Labarsouque and J. Leandri, Nucl. 
Phys. A 613, 342 (1997). 

S. Lemaire, J. Labarsouque and B. Silvestre-Brac, Nucl. 
Phys. A 696, 497 (2001). 

S. Lemaire, J. Labarsouque and B. Silvestre-Brac, Nucl. 
Phys. A 700, 330 (2002). 

S. Lemaire, J. Labarsouque and B. Silvestre-Brac, Nucl. 
Phys. A 714, 265 (2003). 

J. Haidenbauer and G. Krein, Phys. Rev. C 68, 052201 



13 



(2003) . 

F. E. Close, An Introduction to quarks and partons (Aca- 
demic Press, New York, 1979). 

E. M. Henley and G. Krein, Phys. Rev. Lett. 62, 2586 
(1989). 

E. M. Henley and G. Krein, Phys. Lett. B 231, 213 
(1989). 

A. P. Szczepaniak and E. S. Swanson, Phys. Rev. Lett. 
87, 072001 (2001). 

A. P. Szczepaniak, E. S. Swanson, Phys. Rev. D65, 
025012 (2002). 

F. J. Llanes-Estrada, S. R. Cotanch, A. P. Szczepaniak, 

E. S. Swanson, Phys. Rev. C70, 035202 (2004). 

P. Bicudo, J. E. Ribeiro and J. Rodrigues, Phys. Rev. C 
52, 2144 (1995). 

J. R. Finger, J. E. Mandula, Nucl. Phys. B199, 168 
(1982). 

A. Amer, A. Le Yaouanc, L. Oliver, O. Pene, J. c. Raynal, 
Phys. Rev. Lett. 50, 87-90 (1983). 

A. Le Yaouanc, L. Oliver, O. Pene, J. C. Raynal, Phys. 
Lett. B134, 249 (1984). 

A. Le Yaouanc, L. Oliver, O. Pene, J. C. Raynal, Phys. 
Rev. D29, 1233 (1984). 

S. L. Adler, A. C. Davis, Nucl. Phys. B244, 469 (1984). 

A. Kocic, Phys. Rev. D 33, 1785 (1986). 

P. J. d. A. Bicudo and J. E. F. T. Ribeiro, Phys. Rev. D 

42, 1611 (1990); 42, 1625 (1990); 42, 1635 (1990). 

P. J. A. Bicudo, G. Krein, J. E. F. Ribeiro and J. E. Vil- 

late, Phys. Rev. D 45 (1992) 1673. 

H. Mishra and S. P. Misra, Phys. Rev. D 48, 5376 (1993). 

A. Mishra, H. Mishra, V. Sheel, S. P. Misra and 

P. K. Panda, Int. J. Mod. Phys. E 5, 93 (1996). 

F. J. Llanes-Estrada and S. R. Cotanch, Phys. Rev. Lett. 
84, 1102 (2000). 

P. J. A. Bicudo, G. Krein, J. E. F. T. Ribeiro, Phys. Rev. 
C 64, 025202 (2001). 

F. J. Llanes-Estrada and S. R. Cotanch, Nucl. Phys. A 
697, 303 (2002). 

N. Ligterink, E. S. Swanson, Phys. Rev. C69, 025204 

(2004) . 

R. F. Wagcnbrunn, L. Y. .Glozman, Phys. Rev. D75, 
036007 (2007). 

A. Voigt, E. -M. Ilgenfritz, M. Muller-Preussker, 

A. Sternbeck, Phys. Rev. D78, 014501 (2008). 

P. Guo, A. P. Szczepaniak, G. Galata, A. Vassallo, 

E. Santopinto, Phys. Rev. D77, 056005 (2008). 

P. Guo, A. P. Szczepaniak, G. Galata, A. Vassallo, 

E. Santopinto, Phys. Rev. D78, 056003 (2008). 

Y. Nakagawa, A. Voigt, E. -M. Ilgenfritz, M. Muller- 



Preussker, A. Nakamura, T. Saito, A. Sternbeck, H. Toki, 

Phys. Rev. D79, 114504 (2009). 
[59] R. Alkofer, M. Kloker, A. Krassnigg, R. F. Wagenbrunn, 

Phys. Rev. Lett. 96, 022001 (2006). 
[60] A. P. Szczepaniak, private communication. 
[61] T. Nguyen, N. A. Souchlas and P. C. Tandy, AIP Conf. 

Proc. 1361, 142 (2011) . 
[62] C. D. Roberts, M. S. Bhagwat, A. Holl and S. V. Wright, 

Eur. Phys. J. ST 140, 53 (2007). 
[63] R. Alkofer, W. Detmold, C. S. Fischer and P. Maris, 

Phys. Rev. D 70, 014014 (2004) . 
[64] V.N. Gribov, Nucl. Phys. B 139, 1 (1978). 
[65] M. Pak and H. Reinhardt, Phys. Lett. B 707, 566 (2012). 
[66] F. Myhrer and J. Wroldsen, Rev. Mod. Phys. 60, 629 

(1988). 

[67] D. Hadjimichef, G. Krein, S. Szpigel and J. S. da Veiga, 

Annals Phys. 268, 105 (1998) . 
[68] D. Hadjimichef, J. Haidenbauer and G. Krein, Phys. Rev. 

C 63, 035204 (2001). 
[69] A. De Rujula, H. Georgi and S. L. Glashow, Phys. Rev. 

D 12, 147 (1975). 
[70] L. McLerran and R. D. Pisarski, Nucl. Phys. A 796, 83 

(2007) . 

[71] L. Y. .Glozman, Phys. Rev. D 79, 037504 (2009) . 
[72] L. Y. .Glozman, Phys. Rev. D 80, 037701 (2009) . 
[73] C. -Y. Wong, E. S. Swanson and T. Barnes, Phys. Rev. 

C 62, 045201 (2000). 
[74] L. I. Schiff, Quantum Mechanics (McGraw-Hill, New 

York, 1968), pp. 384-387. 
[75] J. Vijande, G. Krein, and A. Valcarce, Eur. Phys. J. A 40, 

89 (2009). 

[76] S.K. Adhikari, Variational principles and the numerical 
solution of scattering problems (Wiley, 1998). 

[77] CNS DAC Services ( http://gwdac.phys.gwu.edu ). 

[78] Z. -w. Lin and C. M. Ko, Phys. Rev. C 62, 034903 (2000). 

[79] Z. -w. Lin, T. G. Di and C. M. Ko, Nucl. Phys. A 689, 
965 (2001). 

[80] A. C. Davis and A. M. Matheson, Nucl. Phys. B 246, 
203 (1984). 

[81] R. Alkofer, P. A. Amundsen and K. Langfeld, Z. Phys. 

C 42, 199 (1989). 
[82] P. Guo and A. P. Szczepaniak, Phys. Rev. D 79, 116006 

(2009) . 

[83] P. M. Lo and E. S. Swanson, Phys. Rev. D 81, 034030 

(2010) . 

[84] B. El-Bennich, G. Krein, L. Chang, C. D. Roberts and 
D. J. Wilson, Phys. Rev. D 85, 031502 (2012). 

[85] C.E. Fontoura, J. Haidenbauer and G. Krein, in prepa- 
ration. 



